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Abstract 

We consider mass correction effects on the polar angular distribution of a 
baryon-antibaryon pair created in the chain decay process e~e~^ — > J/ip — > 
BB, generalizing a previous analysis of Carimalo. We show the relevance of 
the features of the baryon distribution amplitudes and estimate the electro- 
magnetic corrections to the QCD results. 
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I. INTRODUCTION 



In the last few years our understanding of exclusive hadronic processes at high trans- 
fer momentum, in the framework of perturbative QCD, has improved (for a comprehensive 
review and further references see, e.g., ). Theoretical models, based essentially on factor- 
ization ideas, have been elaborated and refined. The main ingredients of these models (which 
from now on we shall indicate as PQCD models) can be summarized as follows: The ampli- 
tude for a given exclusive process is obtained by convoluting two well distinct contributions: 
the first coming from the hard scattering among the partonic constituents of the involved 
hadrons and the second from the subsequent, soft processes which lead to hadronization. 
The hard scattering can be described by means of perturbative QCD techniques, represent- 
ing, in first approximation, each participating hadron by its valence constituents, assumed 
coUinear with the parent hadron and among themselves. As for the soft processes, their 
treatment is outside the possibilities of perturbative methods and alternative approaches 
(hke QCD-sum rules methods or Lattice calculations) are required; in practice, they ap- 
pear in the so called hadronic distribution amplitudes {DA) which describe, for each hadron 
(and independently of the particular process under consideration), how its momentum is 
shared among the valence constituents. Although a complete, formal proof of the validity of 
these factorization procedures is still lacking, at least for exclusive processes, there is enough 
theoretical work which supports these models 

As it should be clear, PQCD models acquire full validity only at very high transfer 
momentum; however, it is not yet clear what that means in practice (that is, at what 
scale we expect these models to become reliable). This is a controversial point 0], even if 
very recently new developments seem to justify the applicability of the models also for not 
so high momenta 0. On the other hand, we must not forget that all the cross sections for 
exclusive processes behave, at high Q^, as an inverse power of . This power increases with 
the number of valence constituents involved in the process and makes more difficult, from an 
experimental point of view, to distinguish and measure these increasingly rare events from 
the bulk of the inclusive processes. We may then summarize the situation as follows: on one 
hand, theoretical models are surely under better control at very high but, on the other 
hand, almost all the experimental information presently at our disposal falls in a range of 

which, while not completely out of reach of PQCD techniques, it is not fully recognized 
as an ideal laboratory for perturbative models. In particular, it is not clear the role played 
by higher order corrections (which are reflected in several possible modifications of the 
basic PQCD models). Unfortunately, the implementation of these higher twist contributions 
is quite intricate. It follows from what we said that a complementary, theoretical and 
phenomenological analysis of all the presently available experimental measurements would 
be of great help in clarifying and improving PQCD models. Thus it is very useful and 
important to undertake all the possible efforts in order to shed light on controversial points 
and improve our understanding of exclusive processes. 

We should also bear in mind that presently only for a few, relatively simple processes, 
calculations have been performed. This is due to the increasing complexity of calculations 
when more and more hadrons (and, as a consequence, partonic constituents) are involved in 
the process. When compared with the experimental results (all of which, with the possible 
exception of the proton form factor, are at intermediate values of Q"^), these calculations 
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show several successes but also some failures. Most of these failures can be attributed to 
violations of the so called helicity selection rules [Q (which are a specific property of PQCD 
models, valid to all orders in the strong coupling constant perturbative expansion) and their 
overcoming requires the introduction of higher twist effects. 

Several attempts have been made in order to implement the original PQCD models 
taking into account higher order contributions, like higher valence Fock states |0, transverse 
momentum effects 0, L 7^ angular momentum components in hadron wave functions [0, 
constituent quarks mass effects |]8|-|10|, diquark correlations inside baryons ||11|| . 



For constituent quark mass effects, in particular, a number of calculations p-pi![] have 
been performed for several experimentally observed processes, some of which are allowed in 
PQCD models while others are forbidden by the helicity selection rules. In these calcula- 
tions the elementary hadron constituents are given a mass, rrii = Ximn, where mn is the 
hadron mass and Xj is the (light-cone) fraction of the hadron momentum carried by the i-th 
constituent, opportunely weighted (in the convolution integral) by the corresponding hadron 
distribution amplitude. This effective mass could take into account (in a global way) several 
higher order effects which have been neglected in the ordinary PQCD models. 

Although from a formal point of view this approach requires further justifications, it offers 
a relatively simple, parameter free, means for practical calculations. The results obtained 
this way can be compared with the lowest order PQCD results, and, as we shall also see in 
the following, lead to predictions for several effects that can be experimentally tested at the 



present time or in the near future (see also ref. ||10|| ). 

Few years ago Carimalo [jl2| considered mass corrections effects on the polar angular 
distribution of baryon-antibaryon pairs produced in the exclusive decay of the J/ip. As it 
can easily be seen, PQCD models predict, due to the helicity selection rules, a distribution 
of the type 1 + cos^ 9g |T|], where 9g is the polar angle which specifies the direction of motion 
of the produced baryon in the J/ip rest -frame. Although experimental results are available 
only for a few baryons JT3| , P^ and are in some cases affected by large statistical errors, there 
are clear indications that the angular distributions behave rather as l + a^ cos^ 6^, with, e.g., 



ttp = 0.62 ±0.11, oa = 0.62 ±0.22 |Tj]. Here (< 1) is a factor which can be expressed from 



the helicity amplitudes for the decay process, as we shall see in detail in the next Section. 

As shown in ref. |T^, mass corrections can in principle explain why < 1 and lead to 
a better agreement between theoretical predictions and experimental results. Furthermore, 
since the parameter is given as a ratio of squared helicity amplitudes, it is independent 
of the exact value of the baryon decay constant, and of several "fine tuning" details of the 
models (for example, how to treat the strong coupling constant in the convolution integrals); 
as such, their effects may be neglected almost completely, as will be clarified by our explicit 
calculations. 

The results of ref. ||12| are obtained using a non-relativistic bound-state model both 
for the decaying J/ ip and the produced baryons. However, while this approximation is well 
grounded for heavy quark bound states, like the J/ip,it can be questionable for light hadrons. 
In ref. [|r^] the proton decay constant was consistently fixed in order to reproduce the J/tp 
pp decay width. A different, widely used approach (see e.g. ref. flp^Sj and references therein) 
is to tentatively fix once and for all the baryon decay constant from QCD sum rules and 
consider different modelizations for the baryon distribution amplitudes, including the non- 
relativistic, QCD sum rules inspired and asymptotic ones. A lattice calculation of the baryon 
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decay constant in the nucleon case is consistent, within the inherent systematic uncertainties 



of the models, with that of QCD sum rules []16|. In this context, the non-relativistic DA 



seems to systematically underestimate absolute quantities {i.e., not obtained as ratios of 
amplitudes) like the decay widths for charmonium exclusive decays, or the hadron form 
factors, by two or even three orders of magnitude, when compared to available experimental 
results (see, e.g., ref. It is then quite reasonable to expect that the use of more refined 

distribution amplitudes could lead to significant modifications in the predicted values for 

There is also a more specific reason to believe that the dependence of on the DA^s may 
be not negligible: the expression of is given as a function of squared helicity amplitudes 
with different values of the constituent helicities. Then, the use of distribution amplitudes 
which, like those derived from QCD sum rules, seem to indicate an unusual sharing of 
the hadron momentum among its constituents, can substantially modify different helicity 
amplitudes. 

Based on these motivations, in the rest of this paper we shall present a new derivation 
of the parameter a^, generalizing the results of ref. to the case of a generic hadron 
distribution amplitude. The only restriction on the DA is that it must satisfy general 
symmetry properties. This allows us to compare the results obtained using different DA^s 
and possibly to get useful informations on them, independently of the way the baryon decay 
constant has been fixed. 

We wish also to recall that an alternative approach for the calculation of has been 
proposed by Parisi and Kada. In ref. Jl^ these authors evaluated the strong contribution to 
ttg for the octet baryons, in the framework of a quark- ( scalar )diquark model for the baryon 
structure. Diquarks are another possible way of taking into account higher order corrections 
(in particular correlations between two valence quarks in the baryon). They have been 
applied, with good success, to several exclusive processes at intermediate values of ^1 



It was stressed by Carimalo and by Glashow et al. []TB[ and will also be argued in the 
following, that electromagnetic (e.m.) corrections to might be by no means negligible. 

The main problem with the e.m. corrections is that they involve the e.m. form factors of 
the octet baryons, for which calculations including mass corrections have not yet been done. 
Only for the nucleon, in the non relativistic approximation, there is a theoretical evaluation 
||19|| . In the other cases we are forced to give for these corrections estimates based on the 
available experimental information. 

This paper is organized as follows: in Section II we derive a general expression for the 
angular distribution of the baryons produced in J/ip BB decays. In Section III we 
concentrate on the strong contribution to a^, which controls the BB angular distribution, 
giving a derivation of the helicity amplitudes required for its calculation. In Section IV we 
discuss in detail the results obtained and their dependence from some subtleties of the models 
that in general, while not modifying qualitatively the conclusions, can substantially affect 
the numerical results, in particular for absolute quantities like decay widths. In Section V 
we give a detailed analysis of electromagnetic corrections, trying to estimate, when possible, 
upper bounds on the consequent overall modification of both from experimental and 
theoretical information. Finally, our conclusions and future perspectives are discussed in 
Section VI. 
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II. DERIVATION OF THE BARYON ANGULAR DISTRIBUTION FOR 

J/t/j BB DECAYS 

Let us consider a. J/ip particle, produced in e~e"'" colliders with unpolarized beams, which 
subsequently decays into a baryon-antibaryon pair: 

e-e+ J/%1) BB . (1) 

The spin density matrix of the J/V", in its rest frame (which is also the cm frame of 
colliding beams, with the electron moving along the positive z direction), has the following 
expression (see, e.g., ref. [pO| ): 

Pmm'(W) = ^ E KnK-,K,K:r,X^-^^, > (2) 
A _,A + 

where the A'^'s are the helicity amplitudes for the process e~e+ J/tp, M is the z component 
of the J/ ijj total angular momentum in its rest frame and Ag- , Xe+ are the helicities of the 
electron and the positron, respectively; is a normalization factor, such that Tr[ p ] = 1. 
It is not difficult to show that: 

p,„,, = if M ^ M' (3) 

where me and are the masses of the electron and the J/V", respectively. As for the second 
step in our process, the decay J/ip BB, we have the general relation pT[| : 



dViJ/i) BB) = — — - 1-4 ^ 



2 \ 1/2 



X E Pmm \^\^\^-M?d^^ , (6) 

where the A's are the helicity amplitudes for the decay of a J/'?/', with third component M 
of the total angular momentum J = 1, into a baryon-antibaryon pair with helicities and 
A^ , respectively; is the mass of the produced baryons. 

Due to the symmetry around the z axis, we can put = in our calculations and 
integrating over yj^ we obtain: 

dV^jg, -> BB) ^ ^_ ( _ 

d{cose^) 8{27r)^\ Ml) 

X E Pao/ I^A.A^;MP • (7) 



5 



Apart from overall factors, independent of 6^^ this is the quantity measured by the 
MARKII and DM2 collaborations Hljll. 

We know from first principles [^] that the amplitudes Axj^\-.m have the following general 
structure 

^a,a^;m(^„V^J = A.A^ <A,-A^(^J exp(^M^J , (8) 

where the "reduced" amplitude ^a^a^ is independent of M and the angular variables and 
the d {6g) are the usual rotation matrices. Using the parity properties for the 

(which imply that A^j^ = v4+_, A = ^4++) and for the spin density matrix p{J/%l)) 

{p-i-i = pi,i), Eq. may be rewritten as follows: 

dTjJ/^ -> BB) ^ _J_ f _ ^rr^\ 
d{cose^) 8(271)^ \ Ml) 

X {|i+_|2(pii + Poo) + 2|i++|Vii} 

X {l + a^cos^^g} , (9) 

where 

{|i+__|2-2|i++|2}(pn-poo) 
a„ = = . (10) 

|A+„|2(pn+poo)+2|A++|2pii 

If, with good approximation (see Eqs. (^,(|^)), we take pn = 1/2 and poo = 0, we finally 
get the simplified expression 



dT{J/iP^BB) 1 f, m 



2 \ 1/2 



rf(cos^^) 16(277)4 \ M^J 

X {|i+^|2 + 2|i++p} {l + a^cos'^^} , (11) 

where 

14 |2 _ 914 |2 
a„ = = . (12) 

We stress again that the results of Eqs. ([TT| ) ,([T^ ) only require the assumptions that the 
one-virtual photon interaction dominates the J/ip production process and that ml/Ml = 0. 

Two further remarks are appropriate: i) As it is clear from Eq. and from the properties 
of the d , 

i+_ = A+_;i(^, =y.^=0) , (13) 

and 

i++ = A++;o(^^, = = 0) . (14) 
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Then it is sufficient to calculate the amplitudes ^++;o in the particular, convenient 

kinematic configuration = V^s = 0, in order to know, with the help of Eq. and of 
the parity symmetry properties, all the amplitudes, ii) As it will be explicitly shown in the 
next Section, if = then A^^ = also; in this case we recover the old PQCD result of 
Brodsky and Lepage that is = 1. 



From Eq. (|g) we easily get the expression of the total decay width for the process: 

2 \ 1/2 

. (15) 



— 1 / TTl 

^{J/^p^BB) = ——{l-4- - 



III. THE STRONG CONTRIBUTION TO THE PARAMETER a^: EVALUATION 

In this Section and in the next one we shall neglect the electromagnetic corrections to 
a^, which will be considered in Section 0, and concentrate on the strong contribution, 
that from now on will be called a^. As we briefiy sketched in Section I, the calculation 
of the helicity amplitudes for the physical process J/V' — > BB, A\ ^ (the suffix s is 

B B ' 

a reminder that we are considering only the strong interaction in what follows), consists 
of several steps. First of all we need to calculate the amplitude for the hard interaction 
among the elementary (valence) constituents of the involved hadrons. To lowest order in 
the strong coupling constant, the only (topologically distinct) Feynman graph is shown in 
Fig. |l| (where the notation is also defined). All other possible graphs of the same order 
can be obtained from this one by a permutation of the final fermionic lines. We do not 
take into account explicitly all these graphs because their contribution is accounted for by 
opportunely choosing the final hadron wave functions. Without giving inessential details of 
the intermediate steps of the calculation, we present directly the expression of this amplitude 
in the particular kinematic configuration Ob = = ^ (the relative momentum between the 
c and c quarks, fc, has also been set equal to zero. This procedure is proper when considering 
L = bound states; for L 7^ the limit should be taken in a subsequent step, i.e. 
after the integration over the angular part of the charmonium wave function). 



A,iA,2A,3,A,-,A,-2A,-3;A.A.(fc = 0; = (^B = 0) = -lQc^g\ 
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Fits ' 



^1 nil \xiyi + {x, - y.fel 



1 



1 



X 



2x11/1 -xi~yi + 2{xi - yifel 2x3^/3 - X3 - 7/3 + 2(^3 - y^fe] 



{ 



xm + X3yi + 2{xi - yi){x3 - y^)^ 



S\g, ^A,2 ,-A,-2 <^A,3 ,-A,-3 Sxg, ,-A,2 ^A,^ ,-A,3 ,xjx,, 



^fl [ "A,i ,Xg^ <JXg^ -A52 "A,3 -A,-3 WA,2 ,A, WA,3 -Ac"Ac,Ac 



+ SXg^ ,Xg^ 5X^2 -Xg2 ^ Xg,^ ,Xg,^ ^A^j .A^f^A^ "A^ 



f^Aqj -A,-, 5a,2 -A,-2 ^Xg,^ ,Aq3 5a,i ,xJXg2 -xJXcXc 

,2 

^B 



X 



X1X3 + yiys - 2{xi - yi){x3 - y3)e\ 

^Xg^ ,^Xg^ 5a,2 ,A,-2 SXg^ - Xg,^ ^ Xg^ , xj Xg^ , - X^^c M + ^fl ^ ^2 '^^2 '^A<73 -A,-3 ^ Xg,^ .A.^A^ ^A^ 



B 

+ ^A^, -Xg^ 5Xg.2 .A92 '^A53 ,A,-3 5Xg^ ,xjXc,-X. 



^B "^Aqj ,Aqj ^Xg2 ,Xg2 ^Xg^ ,Xg^ ^Xc,Xc 



} 



(16) 
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where is the color factor which, once the convolution with the final hadron wave functions 
is made, takes the value = 5/(18a/3) and, as usual, gs = A/47ra7; Xi{yi) represents the 
(light-cone) fraction of the baryon (antibaryon) four-momentum carried by the i-th quark 
(antiquark); = nig/M^. It is clear from this equation the role played by mass corrections 
in allowing spin-flips along the final fermionic lines. We see that higher powers correspond 
to terms where more spin-flips are present. On the contrary, if we neglect quark masses 
Eq. (^) reduces, apart from an inessential constant factor due to a different notation, to 



the result of Brodsky and Lepage |22 



T^°^' = -16c ^ 



We would like to stress that, in principle, the g1 factors coming from different virtual 
gluons should be evaluated at the corresponding values of transferred 0. For the moment 
we do not discuss explicitly this problem, which will be analyzed in detail in the next Section, 
where numerical results are presented. 

Next we convolute the elementary helicity amplitude (for brevity we indicate by {A} 
the collection of all the helicities from which an amplitude depends) with the flnal hadron 
wave functions: 



MI^A^;A.A. = / [di][dy]^B,X^i)Tlx}{i.m„X,{y) , (18) 



where /[rfz] stays for Q dzidz2dz35{l — zi — Z2 — z^) and we used z as a shorthand for 
{zi,Z2,Z3). The amplitude M^^^ refers to the decay of a free c, c quark pair into the flnal 
baryon-antibaryon pair. 

Once the amplitudes have been evaluated (we will do this below), the flnal step 
consists in integrating these amplitudes over the proper cc-bound state wave function, taken 
as usual in the non-relativistic approximation, such that the physical amplitude A\ a- m? 
to which we are interested in, is given by the following general expression: 



/2L + l\^/2 _i 1 



^Ag,A^;A/ - 2^ i ) '-'Ae,-A,A'-^( 



AcA 



S 



iLSJ 



\ An ) ^Ac,-Ac A "-^OAA 



xJd'kMi^^,^.^,^,^ik)Di\i(3,a,0)Mk) , (19) 

where the C's are the Clebsh-Gordan coefficients, A = Ac — Ag, = {k,a,(3) is the relative 
momentum between the c and c quarks and finally ipc{k) is the (momentum-space) charmo- 
nium wave function. In particular, for the J/ip L = 0, so we can take from the beginning 
(in the full non-relativistic approximation) k without loss of generality (this property 
has been used in the derivation of Eq. (p^). Then we can see that: 

A;_;1 = v^7r|i?,(0)|M;_.+_ , (20) 
A'^^.^, = 2n\RM\MU;-,+ , (21) 
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where Rs{0) is the value of the L = charmonium wave function at the origin. By comparing 
the theoretical prediction for the decay width T{J/ip e'e^), V^^ = (16/9)a2|^,(0)|VM|, 
with the available experimental data [^] we can estimate |-Rs(0)| — 0.737 GeV^/^. Here we 
have considered only the two amplitudes that, as it was discussed in the previous Section, 
are sufficient to recover the expressions of all the others, when use is made of Eq. (Rl). In 



Eq. (pri) the relation M^_^_. = M^^._^_, has also been used (the validity of this relation can 



be proved from Eqs. (p!6D,(p!8|)). 

As it is clear from Eq. (plSf ), while the amplitude is the same for all the baryon pairs 
considered, the amplitudes are different for the different baryons, essentially because 
in general the spin-flavor component of their wave function changes. So, we cannot give a 
general, explicit expression for the amplitude M*, but we must consider separately all the 
different cases. The total hadron wave function consists of a color part (which is the same for 
all the baryons and has been englobed in the definition of the color factor of the amplitude 
T''), a spin-fiavor component and a dynamical part which describes, in momentum space, 
how the baryon four-momentum is shared among its valence constituents: the distribution 
amplitude. In general the DA, which as we said previously is a non-perturbative quantity, 
allows for non S'f//(3)-symmetric configurations, as it has been shown by several studies 
performed with the help of QCD sum-rules techniques ||15| , p^ -p^. However, in the particular 



case of an S'f//(3)-symmetric DA, as we shall see, we obtain the same result for all the baryons 
and the only differences are the (experimental) values of the baryon masses . 



In order to explicitly calculate the M amplitudes and their dependence from the DA^s 
we need to take from the literature the available proposed models. QCD sum-rules results 
for octet baryons DA^s exist at present only for the nucleon [p!^ , p^ -|28[| and for the S"*", S~ 



and A |T^. Then, below we give explicitly the expression for the distribution amplitudes 



and for the required amplitudes in these cases. 

The most general wave functions for the baryons considered here are the following: 
for p, n, baryons: 

^B,A,(i) = 2Ab^{¥^b(123)/i,,,(1)/2,-a,(2)/3,a,(3) + y.i,(213)/i,_,,(l)/2,A,(2)/3,A,(3) 

- 2i?BTB(123)/i,A,(l)/2,A,(2)/3,-A,(3) + (1 ^ 3) + (2 ^ 3) } ; (22) 
for the A baryon: 

V'a.aJx) = 2AA^{^A(123)«A.(l)rf-A.(2).A.(3) -^A(213)n_A,(l)rfA.(2)sA,(3) 

- 2RATA{123)ux^il)dx^{2)s-xA^) + all the permutations of (1, 2, 3) } . (23) 
We have introduced the notation 

Rb = ^ , (24) 

where Fb and are constants related to the value of the baryon wave function at the 
origin. In Eq. ( p2|) /i.2,3 are the flavors appropriate to the particular baryon considered 
(/i,2,3 = uud for the proton, udd for the neutron, uus for the S+ and ssd for the S~). 
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Isospin symmetry properties impose several relations between ^g{x) and Tg(a;); in the case 
of the nucleon the relations -R^v = 1, 2TAr(l,2,3) = v9Ar(l,3,2) + v97v(2,3, 1) also hold (see 



ref. |]T5[ for further details). By insertion of these expressions and Eq. (|T6|) in Eq. 
can, after some algebra, derive the following amplitudes. 
For N , baryons: 

+ = / [dx][dy]CB\2DB[pB{l2^)VB{2l^) - 2i?B(^B(321)TB(321) 

+ (/.B(312)rB(132))] +Eb[</^b(132)(/.b(312) -4i?BTB(123)</.B(213)" 
+ 2EBel [v^|(123) + </^|(213) + ^|(312) + 2i?|(2r|(123) + r|(132)) 



we 



Ml 



(25) 



Ml 



__l j [dx][dy]CB[DByB{m + 2R'Bn{m 

+ L'se^[(/5B(132)(^B(312) -4i?sTB(123)¥PB(213)] +2EBej[</^B(123)(/.B(213) 
- Rb{2^b{12?>)Tb{12?,) + ^b(312)Tb(132) + ^B(132)rB(132)) 



For the A baryon 

M', 



/ [dx][dy]CK{ -2Da ^a(123)^a(213) + 2i?A(^A(321)TA(321) 



+ ^a(312)Ta(312))J -EA[</^A(132)y.A(312)+4/?A¥'A(213)TA(213) 
+ 2EKel f</^i(123) + ¥^i(213) + (^1(312) + 27?^ (2X^(123) + T1(312)) 



(26) 



(27) 



Ml 



__A j [di][dy]C^\^D ^[vl{l2?>) + 2RlTl{l?>2) 

- D^tl [(/.A(312)y.A(132) -4i?A</'A(213)TA(123)^ 

- 2EAe; U (123)</.a(213) + 2/?A(</^A(321)rA(321) + y.A(312)rA(312))l \ . (28) 



In Eqs. (]25|)-(p8D the following concise notation has been used: in each if{i, j, k)(p{l, m, n), 
k)T{l,m,n) and ip'^{i,j, k)=ip{i,j, k)(p{i,j, k) product the first term is a function of x, 
the second of y. Furthermore we have defined: 

Db = xiys + xsyi + 2{xi - yi)(x3 - ?/3)e^ 



Eb = 
Cb 



-{xixs + ym) + 2{xi - yi){x3 - y^)t 
2hmTi^al 1 

" ~~v3A^n3?^^rr(^~^^^ 
1 



(29) 
(30) 



X 



X 



2^2 1 



[2xiyi -xi-yi + 2{xi - yi^e 
1 

[2x3?/3 -X3-y3 + 2{x3 - ys^el] 



(31) 
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Let us finally stress that, as it was anticipated, all the M^_^_.^_^_ amplitudes vanish if we 
take rrig (this means that a'^ =1), as it must be, because in this case there is not in the 
model any mechanism which allows for spin flips in the quark-gluon vertices, thus forcing 
the final baryons to have opposite helicities. 

IV. THE STRONG CONTRIBUTION TO THE PARAMETER a^: NUMERICAL 

RESULTS 

In the expressions presented in the previous Section for the amplitudes, from which 
can be evaluated (see Eqs. (p!2D,(p!6D,(p!8D,(|20|), (|2T|) ) we are left with the unknown form of 
the distribution amplitudes, which contain all the hadronization dynamics of the partons of 
the dominant, valence light-cone Fock state. As we said previously, these DA^s are highly 
non-perturbative in nature, so perturbative QCD cannot say about them much more than 
their general, formal solution (expressed as an infinite sum of Appell polynomials with un- 
known, non-perturbative, coefficients) and their evolution in [Q. This is by no means a 
trivial result, but it is not enough to give us reliable expressions for the DA at values of 
presently accessible. For —>■ oo, only the first term in the formal expansion of the wave 
function survives, so we obtain a simple expression, the so called asymptotic distribution am- 
plitude. The asymptotic DA again contains an unknown coefficient, which however may be 
extracted from the best experimental data at hand. The present status of non-perturbative 
methods (namely QCD sum rules and Lattice calculations, in our specific case) allows us 
to calculate only a few low moments of the DA^s. From these, a model expression for the 
DA is proposed by opportunely truncating the infinite, formal expansion of perturbative 
QCD and finding approximate values of the unknown expansion coefficients by fitting the 
calculated moments. There are, of course, several difficulties in this procedure. These are 
mainly refiected in an apparently strong dependence of the model DA from the number of 
available moments; in sizeable differences among the DA^s proposed by different groups; in 
discrepancies between the QCD sum rules and the Lattice results. In our calculation, we 
will take a phenomenological approach, considering all the DA expressions available and 
analyzing how the results for a** depend on each of them. 

Let us first consider the case of the nucleon. There are several versions of QCD sum 
rules model DA^s in the literature. For completeness we give these below, together with the 
non-relativistic and the asymptotic DA's. 




(32) 



ip (x) = 120x1X2X3 . 



(33) 



The DA proposed by Chernyak and Zithnitsky ||2^ 



+ 8.82x1 - 1.68x3 - 2.94] . 
The DA of Chernyak, Ogloblin and Zhitnitsky |T5[ 



(34) 



11 



+ 6.174x2 + 5.88x3 - 7.098]. (35) 



The DA of King and Sachrajda |26 



ip§^{£) = ip'^^{x)[20.16xl + 15.12x^ + 22.68x^ 

- 6.72x3 + 1.68(xi - xa) - 5.04] . (36) 



That of Gari and Stefanis 27 



^%^{x) = (^^^(£)[-1.027x? + 12.307x^ + 25.88x2 

+ 111.32x1X3 + 9.105(xi - X3) - 19.84] . (37) 

The improved (Heterotic) version of the Gari-Stefanis DA, proposed by Stefanis and 
Bergmann [p^ 



!f^^^{x) = (/?^^(x) [-19.773 + 32.756(xi - X3) 
+ 26.569x2 + 16.625x1X3 

- 2.916X? + 75.25x^] . (38) 

In all these cases the nucleon decay constant has roughly the same value: 

= 5.0- lO^^GeV^ (39) 

Furthermore, as we said previously, Rn = 1 and 2T/v(l, 2, 3) = v5Ar(l, 3, 2) + (PnC^, 3, 1) 
(see Eq. (|2|)). 

The only parameter that must be fixed is the value of the strong coupling constant a^. 
For the moment we will limit ourselves to take a fixed value of for all three virtual gluons, 
using the perturbative expression for the running coupling constant: 

"^^^ ^= (lln,-2n;)log(gVA2) ' 

where nc = 3 is the number of colors and ^/ is the number of active flavors (in this context, 
Uf = 4). From this point of view must be set to an overall effective value, pertinent to 
the process considered, and we take the square of the J/ip mass as this effective scale; we 
also use A ^ 0.2 GeV. 

This is by no means the only possible way to treat as we shall see in the following. 
Of course, these ambiguities reflect some uncompletness of the theoretical models. 

In order to evaluate the helicity amplitudes AxgXg-M we finally need to compute quadru- 
ple integrals, which is done numerically. We have tested both our analytical results (the 
expressions of the amplitudes shown in the previous Section) and our numerical calcula- 
tions. For example, it is not difficult to see that if we consider a non-relativistic DA we 
recover, through the different steps of our calculation, the analytical form of proposed 
by Carimalo [|T2|. As for the numerical integrations, we have explicitly checked that if one 
of the usual representations of the Dirac function is inserted in the integrals, and a limit 
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procedure is made for the parameter which enter the representation, the results smoothly 
tend to the analytical ones in the non-relativistic case. 

In Table I we give our results for and for the total decay width V^{J/ip NN) = 
for the nucleon. Proton and neutron results differ in this case only for the small mass 
difference, so we do not present separate results (the situation will of course be different 
when we shall consider, in the next Section, the electromagnetic corrections). 

A comparison of our numerical results for the total decay widths with those of ref. |[T2|JT3 



needs some care. Ref. presents results only in the massless case and using ~ 0.3 and 
|i?g(0)| ~ 0.690 GeV^/^. When opportunely rescaled to these parameterizations our massless 
results are in good agreement with those of ref. |T5[. Note, however, that from an analytical 
point of view our expression for T{J/ip BB) in the massless case differs from that of 
ref. by an overall factor of 80/81. 

In the case of the non relativistic DA our analytical expression of T{J /ip — > BB) agrees 
with that of ref. |]T^. Let us stress, however, that our = 2^/^F/v(ref. |T^); that is, 
F;v(ref. ill) ~ 5 ■ 10"=^ GeV^ corresponds to our F^v ~ 28 ■ lO^^ GeV^. 

As it can be seen from Table I, QCD sum rules DA^s give results which are sizably differ- 
ent from those of the non-relativistic case, even for an observable like which is given as a 
ratio among squared helicities amplitudes. All the results (with the possible exception of the 
GS DA) are in agreement with the experimental measurements, but we must not forget that 
the experimental errors are large at present, and that electromagnetic contributions can in 
principle modify the theoretical values. We see also that there are non negligible differences 
among the various QCD sum rules DA^s. Of course the discrepancies are dramatically larger 
in the case of decay widths, so one can ask why study the (relatively) little dependence of 
on the DA. The point is that, as we shall see better in the following, the same reason which 
makes variations in so small also tend to make freer from several ambiguities of the 
model, giving, in our opinion, a more sensible test for the distribution amplitudes. However, 
in order to discriminate among different DA^s we need better experimental results than those 
presently at disposal. We hope they will be available in the near future. Note also that, 
with F/v given by Eq. (^), QCD sum-rules DA^s give results for the decay widths which are 
in rough agreement with the experiment, unlike the asymptotic and non-relativistic ones 
(here, as in the following, no particular attempt in modifying parameters in order to better 
reproduce experimental results has been made). 

Apart from the nucleon, little is known about the form of the distribution amplitudes of 
the other octet baryons. To our knowledge, the only available models are those proposed by 
Chernyak et al. [|15| for the and A. 



From ref. [T^ we have, for the E 



(^go^(x) = 42(^^^(x)[0.36x^ + 0.24x 



2 



+ 0.14^3 - 0.54xiX2 



-0.16a;3(a;i + a;2)+0.05(xi-a;2)] , (41) 

r^o^(5) = 42^^^(x)[0.32(x? + xl) + 0.16x^ 

- 0.47xiX2 - 0.24x3(xi + X2)] . (42) 

The decay constants are ~ 5.1 ■ 10"^ GeV^, ~ 4.9 ■ lO^^ GeV^. 
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For the S : 



ipi^^{x) = 42(^^^(£)[0.38x? + 0.20x: 



2 



+ 0.16x^-0.26x1X2 



-0.30x3(xi + X2)+0.02(xi-X2)] , (43) 



Tf ^^(x) = 42¥P^^(5)[0.28(x? + xj) + 0.18x^ 

- 0.16xiX2 - 0.35x3(xi + X2)] , (44) 



with |Fh| ~ 5.3 • lO^^ GeV^, \Fl\ ~ 5.4 ■ lO^^ GeV^. 
Finally, for the A: 



ip^^^{x) = 42^^'^(x)[0.44x? + O.OSx^ 
+ 0.34x^ - 0.56x1X2 

- 0.24x3(xi +X2) -0.10(xi -X2)] , (45) 



Tf o^(x) = 42^^^ (x) [1.2{xl - xl) + 1.4(xi - X2)] , (46) 

and |Fa| ^ 6.3 ■ 10-^ GeV^, ~ 6.3 ■ 10"^ GeV^. 

In Table II, we compare the values of a'^ obtained from these QCD sum rules DA^s with 
those obtained in the non-relativistic case and with experimental results (note that, in the 
case of the S"*", experimental data are available for S° only; e.m. corrections can in principle 
be different for E"*" and so a strict comparison of our results with the experimental 
data could be at this stage misleading). We see from Table II that for the S"*", S~ and A 
differences between non-relativistic and QCD sum-rules DA^s results are also more marked. 
This is not unexpected, in that the presence of one or more valence s quarks breaks more 
severely the SUf{3) symmetry implicit in the NR DA. Again, more precise measurements 
of ttg could be very useful in improving our understanding of (at least) the main features 
of the DA^s. By the way, we stress that the decay widths for S^, A seem to be poorly 
reproduced, also when QCD-sum rules DA^s are used. Even if e.m. corrections are not 
accounted for, it is unlikely that they can be able to improve very much the results in this 
sense. 

We want now to analyze better how our results depend on some ingredients of the models 
that, due to their (to some degree) effective nature, are not unambiguously fixed. 

We basically concentrate on the behavior of the results when different ways of treating the 
strong coupling constant are considered. We stress that this is probably (together with the 
value of the baryon decay constant Fb) the main source of indeterminacy of the quantitative 
results of the model, once the DA has been fixed. 

As we anticipated previously, when the calculations of the results quoted in Table I and 
Table II were discussed, the definition of as is in some sense ambiguous. We can make the 
following (reasonable but not derived from first principles) choices: 

i) Consider an overall, effective value of Q^, Ql, which sets the scale to which as must be 
evaluated, using Eq. (^), neglecting that in the model the momenta carried by the 



14 



three virtual gluons in the hard scattering depend on x, y and are different among 
them |l|]. In particular, we choose Ql ~ M^. Slightly different choices have been used 
in the literature (see e.g. ref. this is equivalent to an overall rescaling of the 

decay widths, while is unchanged. 

m) Take a running coupling constant depending on the squared momenta carried by the 
hard gluons, = Q^{x,y) (where i runs over the three gluons) However, if 
this choice is made, we face a serious problem: over the full field of variation of x, y, 
the approach zero. There, then, it does not make sense to speak of perturbative 
calculations, since as increases more and more, invalidating any perturbative power 
expansion. Often people prevents these problems by introducing an ad hoc cut-off, 
Ql, such that 



if > Ql 



^^^^.^^ (lln.-2n,)log(gVA^) (^^^ 



ao = as{Ql) if g2<g, 



where usually one takes = 0.3 or 0.5. In some sense, this is a naive way of accounting 
for non-perturbative effects which prevent as to take higher and higher values (like 
Eq. (pop should imply), even if probably the exact behavior of as is only poorly mimed. 
Even if the — > region is suppressed by the behavior of the distribution amplitudes, 
this procedure is far from being satisfactory. Furthermore it introduces a dependence 
from the new parameter oq. 

in) Take into account in a more serious way non-perturbative effects, which may give rise 
to an effective mass for the gluons (see, e.g., ref. [^). This in turn modifies the gluon 
propagators and the perturbative expression for as, Eq. (|^). To be rigorous, there are 
probably other minor modifications to be introduced into the fermion propagators and 
the quark-gluon vertices as well, if one consistently applies these methods. However, 
since we are only trying to estimate the dependence of our results from different ways 
of treating as, we limit ourselves to consider only the effective gluon mass effects on 
it. In this context, the expression of to be used is of the type: 

a.(Q') = r ' (48) 

(line - 2n/) log (Q2 + 4m2)/A2^ 



where rUg ~ 0.5 GeV ^9 



In Table HI we show, for the nucleon case only (results for the other cases are very 
similar), how the three previous different choices of as influence our results for and the 
decay widths. It is easy to see that is very stable against these changes (as opposed to 
the total decay widths), giving us more confidence on the reliability of our estimates for 
and on their usefulness. 
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Note that the non-relativistic DA has not been taken into account, because in this 
case there is no ambiguity, in that the strong couphng constant is consistently fixed to the 
effective value as{M^). 

Let us conclude by stressing that we also must take into account the evolution of the 
DA with Q^. It is usually said that this evolution has little effect on the results, in that 
its logarithmic behaviour is masked by the stronger, powerlike behaviour (in Q^) of the 
elementary scattering amplitudes. Even if this may seem quite reasonable, we have explicitly 
checked for this assumption, particularly when modified running expressions for have 
been adopted. As a matter of fact, while is almost independent of the DA evolution, 
the modifications of the decay widths due to the evolution are not numerically negligible; 
however, they are not important from a qualitative point of view and do not change the 
results by more than a factor ~ 2, which is within the overall uncertainty we can expect, 
also in the more optimistic hypothesis, for absolute quantities. 



V. ELECTROMAGNETIC CONTRIBUTIONS 

TO ttg 

Up to now we have neglected the electromagnetic contributions to a^. However, it has 



been shown that the corrections due to e.m. processes can be sizable |T2|,|T8[, and a careful 
analysis of cannot avoid considering them. Unfortunately, as it was discussed by Carimalo 



|T^, we are not able, at present, to give a full treatment of these effects. The reason is that 



the time-like e.m. form factors of the baryons are involved, and there is no calculation 
of them including the mass corrections, as it has been done in this paper for the strong 
contribution. Only in the nucleon case, with a non-relativistic DA, this calculation has 



been performed by Ji et al. |T9 



In our treatment we will follow, expand and update the analysis made by Carimalo in 



ref . |]12[ . After a general discussion of these corrections, we shall consider, from two different 
points of view, their estimation (strongly dependent on the baryon e.m. form factors) and 
their contribution to a^, limiting ourselves to the case of the nucleon. In fact, for the other 
octet baryons both experimental and theoretical knowledge of the form factors is too poor 
to allow a sensible phenomenological study. 

We first try to give a full theoretical prediction for in the framework of our model. 
As we said before, at present this is possible only for the nucleon and with a non-relativistic 
distribution amplitude. However, we must not forget that non-relativistic DAs give results 
far from the experimental ones for the e.m. hadron form factors (which are generally un- 
derestimated and given, in some cases, with the wrong sign [|T|). So, we must take the non 
relativistic case as the only one for which at present our model can give a full estimate of 
a^, together with indications on the weight of e.m. contributions, waiting for more reliable 
calculations of mass correction effects to baryon e.m. form factors in the near future. 

Next we shall consider these effects from a more phenomenological point of view, trying 
to estimate their order of magnitude from the available experimental results on the decay 
widths of the J/ip into BB and e~e+ pairs, respectively. 

The reduced amplitudes defined in Eq. (|^), including e.m. contributions, may be written 
as follows 



16 



i+± = + A"^^^ + i;™^ , (49) 

where both e.m. contributions are of order {a /as) with respect to the strong one, and all 
other higher order contributions in {a /as) are neglected. The term emi corresponds, for the 
hard scattering among the valence quarks, to Feynman graphs like those of Fig. [1|, where 
a gluon is replaced by a photon; there are three topologically distinct graphs in this case, 
because we can replace each gluon by a photon. The contribution em2 represents the decay 
of the J /if) into a single virtual photon which in turn is directly (through the time-like e.m. 
baryon form factors) coupled to the final baryon pair. Fig. ^ shows the Feynman graphs 
pertinent to these terms. 

Let us consider, first of all, the emi contribution. It is not difficult to see that, for the 
elementary, hard scattering amplitude, we have 

4 a 

5 as 

where Cj is the electric charge, in units of the proton charge, of the i-th quark, to which 
the virtual photon is coupled, and is given by Eq. (p^. As we said, there are three 
contributions like this, which sum up to give 

4 a 
o as 

Here Qg =61 + 62 + 63 is the baryon electric charge (in units of the proton charge); we 
have also defined the constant 5^ {5b = for neutral baryons and 5b ~ ±10~^ for Q = ±1 
baryons). The subsequent steps which lead to the physical amplitudes A are the same as 
for the strong contribution, so the same relation of Eq. (^Tj) applies also to them: 

i;-=5,i^^ . (52) 

The case of the 6m2 contribution is more subtle. First of all, given that the virtual 
photon couples directly to the final baryons, there is not an equivalent of the hard elementary 
amplitudes T, as in the preceding cases ( or better, this step of the calculation is hidden in 
the baryon form factors). Then we must start from the decay of two free c, c quarks into 
the final baryon pair, which is described by means of the M amplitudes we defined, for the 



strong contribution, in Sect. |ITl 
It can be seen that 



X 



25x,Xg5x,-x,5xg-x- + Sc^^a^a.^a^a^ 
k 



+^ 



Fi{l-Ael)5x^,xA,x,} , (53) 



where kg is the baryon anomalous magnetic moment and are the well known Dirac and 
Pauli baryon form factors respectively, at = M^. 
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Then, from Eqs. ([T|) , (0) , , (|T|) we obtain 



AT^ = ^'n'a\RM\jfGl , (54) 

i;7 = ^7r2a|i?,(0)|^6,Gi , (55) 

where Gfj = Ff + /c^Fg^ and Gf = Ff + {q^ / Arn?^)kBF2 are the Sachs baryon magnetic 
and electric form factors respectively, again at = M|. 

We are now equipped with a formalism for calculating including the leading electro- 
magnetic corrections. Unfortunately a model consistent evaluation requires the analytical 
expressions for the baryon form factors including mass corrections, which at present are 
available only in the particular case of the non-relativistic DA. 

From the paper of Ji et al. [jl9i we find 



G^,(s) = 54F2 7r2a^-(-3-78e2 + 136eJ) , (56) 
Gl{s) = 5AF'^7c'al^i-ll - 246^ + 48eJ) , (57) 



GUs) = 18F'^7rX\{9 + 140e^ - 272e^) , (58) 

G^(s) = 108F>2«,4(7-10e^) . (59) 
It is easy to check that these results reproduce the corresponding massless results,when 



ep,n (see, e.g., ref. I^S] and references therein). However, we point out that mass 



corrections seem to be greater than expected, in that in some cases the higher order terms 
in ep ^ are larger than the leading ones. We never encountered such a situation in our own 
calculations of mass corrections effects. 

By insertion of these results (at s = M|) and Eq. (^) into Eq. (^91) we can estimate 
the value of and its variations with respect to the calculation made only with the strong 
contribution. The results are reported in Table IV, together with the values of the decay 
widths T{J/ip BB), both for the proton and the neutron. We see that percentual 
variations of with respect to the pure strong contribution case are almost negligible: 
the effects of e.m. contributions seem to be much smaller than those due to the choice of 
different DA^s, see Table I. 

Let us consider now what can be said on the e.m. contributions from a phenomenological 
analysis, based upon the experimental information on some ingredients of the model. 

From Eqs. (|g),(|Up,(^, (§§,(§3p we can see that: 



l-p + a-(l + p) 
a„ = , (60) 



B 
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where we have defined 



1 - a* 



1/2 



(62) 



and 



A 



e?Tt2 



1 + 



(63) 



We see then that in order to estimate we need to know and the ratios G^/G^ 



M 



and ^ has been evaluated in Sections fT|,|IV|, while for the two ratios we must resort 

to experimental information. 

Experimental estimates of x+_ can be obtained as follows: from Eqs. ([T5|),(^),(^), ( |55D 
we can write 



+ 2e: 



r^B 



r'B 




(64) 



where F, 



r(j/^ 



BB), r. 



previously defined. From Eq. 



= T{J/tlj e-e+) 
we can evaluate Xj 



and all other quantities have been 
, once 

m. 



\Gfj\ and the ratio G^/G^ 



M 



have been given. We take the experimental data on F^^ , Fee from ref . 

Let us first consider the proton case. Very recently the time-like proton magnetic form 
factor at ~ 10 GeV^ has been measured \Glj{Ml)\ = 0.026 ± 0.002 (experimental 
results on G^,^ were not available at the time ref. [jl2[ was published. By analytic continua- 



rp 



Ml) 



rP 



-Ml) 



0.012. 



tion, it was assumed that to leading order in as 
This does not seem to be the case pO|). 

So, for the proton the only quantity we require which is not experimentally measured 
at present is the ratio |Gf/Gff|. One possibility is to assume approximate validity of the 
empirical relation G'^^/G% ~ 1//Up. However, there is nothing assuring that this relation can 
be verified in the kinematical regime of interest here. For such reason, we will subsequently 
take this ratio as a free parameter, allowing for reasonable variations of its value around the 
expected, empirical one, in order to study the dependence of from this uncertainty. 

Given that our results depend on the values of Gff and |Gf/Gff| separately, let us 



30|. 



clarify how measurements have been taken from ref. 
measured value of \G\j\ depends on the value of the ratio jC^/G^I. 



It is seen from ref. that the 
This is so because 

there is not enough statistics in order to derive separate expressions of G^j and G^^. So, 



in our phenomenological analysis we proceed as follows: once the ratio G^^/G\^ has been 



fixed to some value, we 



can estimate G\^ 



from the experimental data (see Eq. (4) and 
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Table I of ref. [BO]) and using the same fitting procedure of ref. pO| to get the form factors 



at s = The first two columns of Table V give corrected values of and percentual 
variation with respect to the pure strong contribution, when use is made of the empirical 
relation G^^/G\j ~ We can see that e.m. corrections can be relevant, and we probably 

cannot neglect them in a complete evaluation of . We stress however that these estimates 
depend on the theoretical value of the strong contribution (see Eqs. (|60l-^)), as it can be 
seen from the spread of Aop over the various DAs considered. However, since we cover a 
large range of possibilities for the strong contribution, we think the order of magnitude of 
e.m. corrections can be safely estimated to be of 5-10 %, for the proton. These estimates 
must be taken as upper limits. 

Fig. ^ shows the variation of Op as a function of the ratio G^^/G\.j and for some indicative 
DA"?,. The central value correspond to the possible expected behavior, G^^/G^j ~ 

We see from Fig. |^ that our previous conclusions are not greatly modified; only for big 
variations of the behavior of the G^^/G\j ratio from the empirical one, we can have relevant 
differences. Our conclusion is that the full range of variation of with respect to is of 
the order of 0-15 %. 

In the case of the neutron the experimental information is also poorer than for the proton. 
We limit ourselves to extrapolate well known phenomenological behavior of the neutron form 
factors to the kinematic region ~ M^. From this point of view, we can at first assume 
~ and G^.j/jin — G%i/f^p, at = M|. The last two columns of Table V show the 
corresponding results for and its percentual variation with respect to the pure strong 
contribution a^. We see that this time e.m. corrections induce a decrease of because 
is negative. As for the absolute variation, this is slightly larger than in the proton case, and 
in the range of 9-15 %. Even if we do not report here the results, we have analyzed also for 
the neutron how e.m. corrections change when G^j and G^/G'%j are allowed to vary around 
the empirical values given above. We find again a situation similar to that of the proton, 
possibly with a greater spread of variation. 



VI. CONCLUSIONS 

The study of exclusive processes involving hadrons at high energy scales, in the frame- 
work of perturbative QCD models, has made big improvements in the last years. However, 
several problems remain open and questionable, given that almost all the experimental in- 
formation concerning exclusive processes is at intermediate energy scales. In such situation, 
several higher order corrections may not be negligible and can contribute significantly, at 
least for those processes that are forbidden to lowest order in the PQCD models, but are 
experimentally well established. Unfortunately the implementation of these higher twist 
effects is very intricate, and only very recently significant progress has been made. Among 
these contributions, mass corrections effects for the valence constituent quarks of the light 
hadrons involved are very promising. In fact, if from one hand a more rigorous theoret- 
ical explanation of their exact origin is required and auspicable, on the other hand their 
implementation is relatively easy and free of ambiguities: once their contribution is taken 
into account there are no free parameters to be fixed, and everything goes on without any 
further assumptions. A number of exclusive processes (in particular charmonium decays) 
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have been analyzed in this framework in the last years, and several interesting consequences, 
experimentally testable at present or in the near future, have been proposed to check the 
validity of this model as compared to the PQCD ones or to alternative models including 
higher order effects. In this paper we have analyzed the effects of mass corrections on the 
angular distribution of octet baryon pairs created in J/%1) decays. Even if experimental mea- 
surements of these angular distributions have, for the moment, a low statistics (at least in 
some cases) we hope they can be improved in the near future. 

This was also the subject of an earlier paper by Carimalo. However Carimalo used, as 
a first approach, a simplified, non-relativistic model for the produced baryons. Over time it 
has been shown that the behavior of light hadrons, like the octet baryons, probably demands 
more accurate models of the baryon distribution amplitudes, like those we have considered 
here in a more general context. From a theoretical point of view, there are several reasons 
why the study of these angular distributions is interesting. First of all they are governed by 
a parameter which is given as a ratio of helicity amplitudes. As such, this parameter results 
to be quite independent of some details of QCD models presently at our disposal, details 
which are not fully understood and can modify sizably the numerical results (albeit not 
the qualitative ones). Secondly, as we have shown, this parameter is sufficiently sensitive 
to the precise form of the distribution amplitudes to allow a discrimination among their 
main features, as soon as higher precision measurements will be available. Considering 
all the available models for the distribution amplitudes (mainly based on QCD sum rules 
calculations) we have shown indeed that the spread on the values due to the change of 
the DA is of the order of 10-20 %, in the nucleon case. Little is known about the DA'S of the 
other octet baryons. Whenever possible, our calculations show an even larger dependence 
of from the exact form of the DA, varying in the range of 20-50 %. Therefore a more 
accurate measurement of could allow us to discriminate among different model DA'S. 

We have also considered in detail the role played by electromagnetic corrections. In fact, 
we cannot evaluate exactly these corrections at present. However, we give reasonable upper 
limit estimates of these contributions in the nucleon case, making use of all the possible 
available informations , both theoretical and experimental. 

Theoretical estimates (using a non-relativistic approximation for the final baryon dis- 
tribution amplitudes) suggest that e.m. corrections could be negligible, both in the proton 
and neutron case. 

On the contrary, we estimate from a phenomenological analysis that e.m. corrections 
could be of the order of 5-10 % for the proton and of 10-15% for the neutron. In the proton 
case e.m. corrections tend to increase the value of , while for the neutron they produce the 
opposite effect. In both cases, these corrections could be comparable to the modifications 
induced in the pure strong contribution by different choices of the DA. 

We can say very little about e.m. corrections for the other octet baryons, since little is 
known about them, both experimentally and theoretically. So we limit ourselves to report, 
when possible, the pure strong contribution. 

Of course, only more precise experimental information on the e.m. baryon form factors 
could allow to evaluate the exact contribution of e.m. effects. Alternatively, the form factors 
could be evaluated in the framework of our model (including mass corrections), and a more 
consistent, theoretical result could be obtained. We leave this as a subject for future work. 
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FIGURES 



FIG. 1. The Feynman diagram which, to lowest order in Og, describes the elementary process 
QQ ~^ Q1Q2Q3Q1Q2Q3, for a quarkonium state with charge conjugation C = —1. In the QQ cen- 
ter-of-mass frame, = {E, k/2) and = {E, —k/2), where k is the relative momentum between 
the c and c quarks; qi = xips and qi = yiPj^ {i = 1)2,3), with = {E,pb), p'^ = {E, —ps) and 
Pb = {psm9BCOSipB,psiyi6BsinipB,pcos9B)- a,b,c, i, j, 1,1' ,mi^2,3i ^^1,2,3 are color indices; the 
A's label helicities. 

FIG. 2. The two Feynman graphs which describe the leading (in a/os) electromagnetic 
corrections to the QCD lowest order term (see Fig. ||): a) The contribution emi, obtained by 
substituting one of the virtual gluons in Fig. |^ with a photon. There are two other contributions, 
obtained by replacing each of the other two gluons with a photon (see Fig. || for the notation), b) 
The contribution em2, coming from the direct coupling of a virtual photon produced in the QQ 
decay with the final baryons. 

FIG. 3. The dependence of the full (QCD+QED) angular distribution parameter (see 
Eq. (|l2|)) from the poorly known Ge/Gm ratio, in the proton case. The vertical, solid line corre- 
sponds to the well known empirical behavior G^/G^ ~ 1//Up. Different lines correspond to the 
various distribution amplitudes considered for the strong contribution (see text): NR (solid); ASY 
(dashed); COZ (dotted); KS (dot-dashed); GS (double dot-dashed). In this plot, the GZ and 
HET DA results are almost indistinguishable from the GOZ and NR ones, respectively. 
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TABLES 



TABLE L The strong contribution to and to the decay width T{J/iIj BB) for the 
nucleon. Results obtained using the different nucleon DA^s considered in the text are compared to 
experimental data. 



_DA aj^ r-^^ • lO^(GeV) 

NR 0.688 0.002 

AS 0.667 0.026 

CZ 0.561 0.587 

COZ 0.565 0.826 

KS 0.591 1.255 

GS 0.963 0.168 

HET 0.689 1.671 

MK2^ 0.61 ± 0.23 1.85 ± 0.27 

DM2*^ 0.62 ± 0.11 1.63 ± 0.38 



^MARKII Collaboration, ref. 
^DM2 Collaboration, ref. Q. 

TABLE II. The strong contribution to and to the decay width T{J /ip — > BB) for the S"*", 
and A. Only the non-relativistic and the COZ DA^s are available for these particles. 



H- A 



DA 


a% 


(GeV) 


a% 


■ 10' 

(GeV) 


a% 


■ 10' 

(GeV) 


NR 


0.431 


0.002 


0.274 


0.002 


0.513 


0.003 


AS 


0.417 


0.032 


0.265 


0.016 


0.497 


0.017 


COZ 


0.687 


58.151 


0.537 


45.519 


0.770 


7.252 


MK2'^ 


0.7 ± 1.1 


1.35 ± 0.35 


-0.13 ± 0.55 


0.97 ± 0.25 


0.72 ± 0.36 


1.35 ± 0.27 


DM2'= 


0.22 ± 0.31 


0.91 ± 0.27 




0.60 ± 0.15 


0.62 ± 0.22 


1.18 ± 0.26 



^Experimental results are available only for the J/tp ^ S'^S'^ case. 
i^MARKII Cohaboration, ref. ||]. 
'=DM2 Collaboration, ref. Q. 
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TABLE III. Dependence of the strong contribution to and to T{J /ip — > BB) on different 
behaviors of the strong couphng constant inside the convolution integral of Eq. ([l8|), in the 
proton case. First column is the same as Tab. I, with from Eq. (^) at = M^; in the 2nd and 



3rd column use is made of Eq. (|47D , with oq = 0.3 and 0.5, respectively; the 4th column presents 
the results when Eq. (H^ is used. 





as = 


0.275 




< 0.3 


as 


< 0.5 


as 


(ml) 


DA 




(GeV) 




(GeV) 


a% 


^■nn ■ 10' 
(GeV) 


a% 


T'' - ■ lO'^ 
(GeV) 


AS 


0.667 


0.026 


0.667 


0.044 


0.666 


0.466 


0.667 


0.210 


cz 


0.561 


0.587 


0.561 


1.144 


0.564 


10.013 


0.561 


5.716 


coz 


0.565 


0.826 


0.565 


1.369 


0.567 


12.078 


0.564 


6.894 


KS 


0.591 


1.255 


0.591 


2.076 


0.592 


18.343 


0.591 


10.549 


GS 


0.963 


0.168 


0.963 


0.275 


0.954 


1.917 


0.963 


1.160 


HET 


0.689 


1.671 


0.686 


2.763 


0.682 


24.245 


0.691 


14.020 



TABLE IV. Theoretical predictions for and the decay width T^^(J/ip — > BB) , including 
e.m. corrections, for the proton and the neutron. The 2nd column gives the percentual variation 
of Qg with respect to the pure strong contribution. Only the non-relativistic DA is considered, 
and the results of Ji et al. |l^ for the nucleon e.m. form factors are used. 





IN 


Aa^(%) 


Ljv^v • lO^(GeV) 


p 


0.696 


1.2 


0.002 


n 


0.677 


-1.4 


0.002 



TABLE V. Results for Up^n obtained by adding to the theoretical values of a^^^ phenomenolog- 
ical estimates of the e.m. corrections. All the proposed DA^s are considered. The G^(M^) and 
G^"(M|)/G^j"(M0) values are fixed by using, when available, experimental data or by extrapo- 
lating them to the kinematical region under study (see text for details). 







P 






n 




DA 


Op 




Aap{%) 


a„ 




Aa„(%) 


NR 


0.725 




5.3 


0.628 




-8.9 


AS 


0.708 




5.9 


0.605 




-9.7 


CZ 


0.620 




10.1 


0.483 




-14.7 


COZ 


0.623 




9.9 


0.487 




-14.6 


KS 


0.645 




8.8 


0.517 




-13.2 


GS 


0.956 




-0.7 


0.954 




-0.8 


HET 


0.726 




5.3 


0.629 




-8.9 



26 



P5 



O O O O O 

• ••••• 





PcAcJ 



'91 

mi 



ni 



92 

m2 



*93 



?b-^9i ?2,Ag2 qs^K. 



'93 



Fig. 1 




nil 



Til 



1712 



n2 



93 

m3 



93 



Fig. 2a 




Fig. 2b 



